A dynamical description of neutron star crusts 
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Abstract. Neutron Stars are natural laboratories where fundamental properties of matter 
under extreme conditions can be explored. Modern nuclear physics input as well as many-body 
theories are valuable tools which may allow us to improve our understanding of the physics of 
those compact objects. 

In this work the occurrence of exotic structures in the outermost layers of neutron stars 
is investigated within the framework of a microscopic model. In this approach the nucleonic 
dynamics is described by a time-dependent mean held approach at around zero temperature. 
Starting from an initial crystalline lattice of nuclei at subnuclear densities the system evolves 
toward a manifold of self-organized structures with different shapes and similar energies. These 
structures are studied in terms of a phase diagram in density and the corresponding sensitivity 
to the isospin-dependent part of the equation of state and to the isotopic composition is 
investigated. 



1. Introduction 

Neutron stars (NS) are known to be born after the core collapse of supernovae explosions. After 
its birth the NS has cooled below 10 9 K and a solid crust is supposed to be formed in the 
outermost layers. Even if the crust contains only a small part of the star, it plays an important 
role in its evolution. Indeed, the properties of the crust, notably the transport properties, can 
strongly influence processes like cooling, accretion and X-ray bursts. 

Heavy ions and the crust of NS share basic many-body aspects due to the fact that they are 
made of interacting nucleons. Having more neutrons than protons, many physical properties of 
both kind of objects are related with the asymmetry dependence of the nuclear interaction. 

In this work we were interested in dynamical aspects of NS crusts, in particular, with the 
occurrence of exotic structures. In that region, at densities just below the saturation value, the 
nuclear matter is said to be "frustrated" : subjected to the competition of the short-range nuclear 
attraction and the long-range Coulomb repulsion. As a consequence of this competition nuclei 
may organize forming non-spherical composites. Structures such as rod-like and slab-like nuclei, 
the so called nuclear "pasta", were initially predicted to exist in the ground state of nuclear 
matter at subnuclear densities by static approaches [1]. Assuming well denned symmetries, 



these models suggested the appearance of these structures in a well defined order from spherical 
to cylindrical, planar, complex structures and finally to uniform matter, with increasing density. 
This typical sorting of structures was confirmed by recent works [2] and suggested that pasta 
phases can be formed in the cooling process of hot NS and during the compression stages taking 
place during the collapse of supernova cores [3j. In both cases it is worthwhile to give a dynamical 
description of pasta formation. 

The interest in structural phases in NS crusts is double. On one side, their existence may have 
various astrophysical consequences. Their presence can influence the scattering of neutrinos and 
modify the neutrino trapping with the corresponding consequences on the processes of supernova 
core collapse or NS cooling. Another aspect which can be sensitive to pasta formation are the 
elastic properties of the star (namely those phenomena depending on the crust-core boundary 
properties) which may influence gravitational wave radiation. On the other side, from the nuclear 
physics point of view, their formation can be connected with fragment formation in heavy ion 
collisions. These phenomena share many aspects since they are both driven by volume, Coulomb 
and surface energies. In particular, the sensitivity to the nuclear equation of state has been 
suggested [I] as in cluster formation. 

In this work we present a model describing the dynamics of nuclear matter in the crust of 
a NS. It is an extension of a dynamical approach initially developed in order to describe heavy 
ion collisions at intermediate energies [5j. In contrast to the many previous studies employing 
static frameworks, this model allows to simulate the dynamical processes in inhomogeneous 
nuclear matter using a large number of nucleons without any assumptions on the structure of 
nuclear matter. Since these structures involve a large number of low energy configurations, 
in this work a pure mean field description of the nuclear dynamics has been performed. In 
this approach, starting from initial crystalline lattices of nuclei, with different symmetries, non 
spherical structures occur as the result of microscopic self-organization processes. The survival 
of those meta-stable equilibrium structures has been checked over thousands fm/c [6]. In this 
work, the influence of the equation of state (EOS), of the isotopic compositions and of the lattice 
perturbations on the formation of exotic structures are sudied. 

This contribution is organized as follows. In Section 2 the bases of the model are presented. 
In Section 3 we analyze the occurence of structures in Oxygen lattices with different isotopic 
composition. The effects of the EOS and the response to lattice perturbation are discussed. In 
Section 4 our concluding remarks are given. 

2. A dynamical model for the neutron star crust 

At temperatures T < 1 MeV and densities between 10 6 g/cm 3 at the surface to 3xl0 13 g/cm 3 
at the interface with the core, the stellar matter can be modelled by a neutral mixture of 
nuclei, electrons and free neutrons. Nuclei are expected to form a crystalline lattice immersed 
in a degenerate relativistic gas of electrons, which may be described by a uniform density 
distribution [2j. In order to describe the microscopic dynamics of the nuclear many-body system 
we have developed a dynamical approach, which is based on the DYWAN model [5] for heavy-ion 
collisions at intermediate energies. This approach essentially gives a prescription for solving an 
extended time-dependent Hartree-Fock (ETDHF) equation for the one-body density matrix, of 
the form: 



where p is the one-body density matrix, h(p) is the one-body Hamiltonian and )C(p) is the 
collision term, providing the irreversible evolution of the system due to the coupling of p with 
higher order multiparticle correlations. In the present calculations we are essentially concerned 
with low energy configurations of the system. We will suppose that transitions towards excited 
single particle states are inhibited and that the dynamics is mainly governed by the mean field. 



ihp= [h(p),p] +/C(p). 




Equation (pQ), with IC(p)=0, is solved by projection of the one body space of states onto 
a convenient basis. These functions are splines [8], a particular case of wavelets, which are 
mathematical (non analytical) objects with many interesting properties and which behave like 
a moving basis. 

In a first step, a static Hartree-Fock (HF) self-consistent procedure is implemented in order to 
get nuclear composites either in their ground states or in excited states according to mechanical 
or thermal constraints. With these nuclei, prepared in a single Wigner-Seitz (WS) cell, a 
super-cell is constructed by stacking unitary cells in three dimensions. Their number and 
symmetries are different depending on the kind of lattice we are interested in. The complete 
three-dimensional lattice is built up from the super-cell under periodic boundary conditions. In 
this work simple cubic cells (SCC) have been considered, but other symmetries can be arbitrarily 
chosen. 

For this calculation we have chosen, for simplicity, a density-dependent zero-range effective 
interaction, with the following self-consistent field: 

y q (M) - — p + p + —^t + t-2«+ (2) 
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where p n and p p stand for neutron and proton densities, p = p n + p p , £ = p n — p p , q=l/2 
for neutrons and -1/2 for protons, poo =0.145 fm -3 is the saturation density of infinite nuclear 
matter and V q is the Coulomb potential. The calculation of V q is performed using the Ewald 
summation techniques [S], which is adapted to the calculation of long range potentials in periodic 
systems. It consists basically in recasting the Coulomb potential into two convenient terms, each 
one of them can be calculated in a fast and efficient way [B] . 

The parameters of the force are related with those corresponding to the usual Skyrme force, 
to, ts, xq and 23, by the following relationships: 

3 Poo ., 3(V + 2K+\ 
*° ~ ~1~* t3 ~ 4^ t3 ' 

with xo = X3 = —1/2, i/=l/6, and: 

t' /poc = -356MeV fm 3 t' 3 / ' p u + l = 303MeV fm 3 ^^ 1 ) . 

The parameters c and £1 are related to the coefficients of the mass formula J and L corresponding 
to, respectively, the volume-symmetry and to the density-dependent symmetry energy as follows: 



6M' 

L = 6c + tt + L kin 



Lki n is the kinetic energy contribution to L. These parameters are still uncertain because they 
are not completely constrained from nuclear data, then a study of their influence on the overall 
dynamics seems necessary. In this work c has been fixed to 20 MeV and SI varied in the range 
-125 and -75 MeV in order to analyze the sensitivity of the dynamics to this interaction. The 
density dependence of the symmetry energy has been determined here in a phenomenological 
way according to its current estimate at the saturation value in pure neutron matter. The values 
of the parameters reproduce the principal static characteristics of nuclei, as binding energies, 
radii and equilibrium densities. The associated incompressibility modulus in symmetric matter 
^00= 200 MeV corresponds to a "soft" EOS. 



The energy density per baryon ui is defined as: 
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where i^kin corresponds to the kinetic contribution. In Fig. Q] the calculated values of ui in 
pure neutron matter are represented as a function of the mean density for different values of 
the asymmetry parameter fi. Together with our calculations, the results of Ref. [ID], in the 
framework of a microscopic variational model, are plotted in triangles. 




Figure 1. Neutron energy density per baryon as a function of the density at different values of 
the asymmetry parameter O. The results of Ref. [ID] are in triangles. 

Following the iterative procedure, Oxygen nuclei have been prepared in WS cells with the 
force given by Eq. ([2]) . The three-dimensional lattice is then constructed with replicas of basic 
SCC with periodic boudary conditions. In this case the supercell is composed by 27 WS cells. 
In Fig. [2] is represented the projected density profiles on the symmetry plane (X,Y) of a SCC 
lattice of Oxygen isotopes with proton fraction x=0.5 and mean density (p}=0.02 fm~ 3 . Other 
kind of lattices with different symmetries and species have been considered in Ref. . 

As stated before, once the initial conditions are given, the system is led to evolve according 
to Eq. ([!]) turning off the collision term. The resulting equation is solved by projection of single 
particle wave functions onto a spline wavelet basis a(f,t) [8\. These functions are orthogonal, 
they present defined symmetries and compact support and, even if they are not analytical they 
can be approximated by simple analytical forms. At the same time, they provide strongly 
compressed representations of the system with an accurate description of single particle wave 
functions. 

In configuration space wavelets are functions of a set of correlated generalized coordinates 
{£,,X,tt,4>}, which correspond to first and second moments. In the x dimension they read : 

< X >= < (x- i x ) 2 >= Xx < Px >= K x < (p x - 7T X ) 2 >= 4> x , 

with similar expressions for y and z components. 

It can be shown [6j that wavelets satisfy modified TDHF-like equations, while their 
corresponding occupation numbers evolve through a master equation. The equations of motion 
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Figure 2. Density profiles of ie O lattice on the (X,Y) plane. 



for the generalized coordinates can be obtained from a viariational principle [§]: 
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V being the wavelet transform of the effective nuclear potential V =< a\V \a >. It is 
important to underline that the self-consistent mean-field is a function of the density extended 
to the overall supercell and not only to a single WS cell. 

Equations [3j|6l which solve the TDHF equation for the nucleonic system, are Hamilton-like 
equations for the centroi'ds and the widths of the moving basis of wavelets. 

3. From crystalline to exotic structures 

Let us now consider the initial arrangement of Oxygen isotopes, slightly deformed along the z 
axis, in a SCC lattice with proton fraction x= 0.1 and mean density (p) =0.072 fm -3 . Some 
snapshots of the evolution are shown in Fig. [3l The system goes through some preferred 
structures, in this particular case, it oscillates between rods and slab-like shapes. These plots 
represent points in configuration space for which the density is higher than a reference value: the 
threshold density pt- The corresponding value in Fig. [3]is 0.065 fm" 3 . For a given configuration 
the choice of this quantity determines the underlying structures. In consequence the occurrence 
of a given shape depends not only on the mean density but also on the choice of pt- Even more, 
for a given choice of both quantities, different embedded structures can be found. 

In order to characterize unambiguously these structures we have utilized the principles of 
morphology recognition from Morphological Image Analysis (MIA) techniques, which assign 
numbers to the different shapes according to some geometrical prescription [12j . 
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Figure 3. Time evolution of neutron density in oxygen SCC lattice with x p =0.1, (p}=0.072 
fm~ 3 and threshold density 0.065 fm~ 3 . 

Plotting the different structures in the density plane {pt vs (p)) gives an overall view of 
the most probable structures occurrence for different initial configurations. In Fig [5] the 
corresponding diagrams for Oxygen isotopes SCC lattices is represented in the asymmetric 
x=0.2 case, for two values of the asymmetry parameter $7: -120 MeV (a) and -75 MeV (b). 
The different regions in grey scale correspond to the various structures, from black colour for 
spherical nuclei to white colour for sponge-like shapes. In any case, for a given threshold density, 
the five standard types of pasta phases emerge naturally and in the same order with increasing 
density, as predicted by static models [1] and confirmed by recent calculations [2J. 




Figure 4. Neutron threshold density versus the neutron mean density normalized to the 
saturation value for two values of Q and for proton fractions x=0.2. 



In this example the system experiences some sensitivity to the asymmetry-dependent 
contribution to the effective force. Indeed, in Fig. 2] the stiffer force, corresponding to $7=-75 
MeV, is shown to favour the occurrence of spherical and sponge-like structures, while restraining 
the appearance of cylinder and slabs. This is consistent with the fact that in the stiff case the 
potential is more repulsive, namely at high mean densities. The consequence is that more and 
more neutrons are dripped or localized at the surface. In the last case, quasi-free neutrons 
contribute either to build heavy spherical nuclei or to link neighbouring clusters is all directions, 
giving rise to the so-called sponge-like structures. 



4. Effects of lattice perturbations 

As underlined by other authors |13| the inner crust of a neutron star is an extremely complex 
system in which the appearance of disordered phases can be favored. There is a special interest 
to investigate how nuclear matter, for a given proton fraction, self-organizes dynamically with 
the introduction of slight perturbations of the initial conditions. In order to study the feasibility 
of those disordered phases, we considered SCC lattices of oxygen isotopes, with x p =0.2 and 0.5, 
the positions of which have been initially slightly shifted at random from the lattice sites. In Fig. 
0the snapshots of the supercell density at t=600 fm/c is represented in the cases in which the 
proton fractions are x p =0.5 (a) and 0.2 (b). In both cases the average density is < p >= OAp^, 
the threshold densities pt are 0.04 and 0.05, respectively. 

Starting from similar initial conditions the system evolves towards unlike configurations 



Time=600 fm/c Time=600 fm/c 




Figure 5. Snapshots of neutron density at t=600 fm/c in perturbed oxygen SCC lattice with 
(p}=0.06 fm~ 3 for (a) x p =0.5 and (b) x p =0.2. 

depending on their corresponding proton fraction. In the symmetric case the system organizes 
in several intermediate mass clusters, while in the asymmetric case it evolves towards a unique 
massive cluster occupying the entire supercell. In any case, a gradual loss of the initial 
symmetries is observed in both pictures. From an ordered array of regular deformed nuclei 
the system goes through a heterogeneous configuration of clusters with different masses and 
shapes. It is worthwhile to underline that in Ref. [6] the model has been shown to be extremely 
well conditioned, conserving metastable configurations during time intervals of several thousands 
of fm/c. This fact permitted to check the corresponding numerical accuracy and reject possible 
numerical instabilities. 

The small initial perturbations then generate sufficiently strong density fluctuations so as 
to provoke a different rearrangement of matter molding heavy aggregates. The fact that at 
lower proton fraction a single (infinite) aggregate occurs can be interpreted as a consequence 
of the increase of neutron wave functions diffusivity. This is a direct consequence of wavelet 
spreading which results from the dynamical evolution. Indeed, as the proton fraction decreases, 
the number of dripped neutron increases and nucleons become more delocalized with significant 
overlap in configuration space. The competition between nuclear and Coulomb interactions, 
which is local, covers broader spatial regions in this case. 

5. Conclusions 

In this work we present a non dissipative, self-consistent, dynamical description of nuclear matter 
in conditions of density and temperature similar to those which should exist in the outermost 
layers of neutron stars. 

In this description the small excitation energy which is initially deposited in deformation allows 
the system to explore a landscape of structures. The occurrence of most probable shapes is 



analyzed by means of structure diagrams in the threshold density versus mean density plane. 
Not only the ordered five standard types of pasta emerge naturally but also other intermediate 
shapes between them which can be characterized through morphological analysis techniques. 
Phase diagrams are shown to be sensitive to a variation in the intensity of the isovector part 
of the self-consistent potential. The observed trend is the increase of spherical or sponge-like 
structures with the stiffness of the force. 

In the present pure mean-field description lattice symmetries are preserved along the temporal 
evolution of the system. In all cases, the effects of perturbing the initial lattice is to break 
those symmetries. Nevertheless the characteristics of the final distribution depend on the 
isotopic composition. Symmetric systems organize in several intermediate mass structures, while 
asymmetric ones collapse in a unique massive cluster occupying the entire supercell. 

In this work a simple phenomenological local potential has been implemented. The formation 
and transitions between the observed non-equilibrium structures can be strongly modified using 
more complex non-local effective forces and beyond a pure mean-field description. The study of 
the transport properties of the crust can be then performed in more reliable conditions. These 
investigations are currently in progress. 
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